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Abstract. We apply the similarity renormalization group (SRG) approach to evolve 
a nucleon-nucleon (A'A') interaction in leading-order (LO) chiral effective field theory 
(ChEFT), renormalized within the framework of the subtracted kernel method (SKM). We 
derive a fixed-point interaction and show the renormalization group (RG) invariance in the 
SKM approach. We also compare the evolution of A'A' potentials with the subtraction scale 
through a SKM RG equation in the form of a non-relativistic Callan-Symanzik (NRCS) 
equation and the evolution with the similarity cutoff through the SRG transformation. 

1 Introduction 

The standard method for the non-perturbative renormalization of nucleon-nucleon (A^A^) interactions 
in the context of chiral effective field theory (ChEFT), which is inspired by Wilson's renormalization 
group LU, consists of two steps II2I3II . The first step is to solve the Lippmann-Schwinger (LS) equation 
with the A^A^ potential truncated at a given order in the ChEFT expansion, which consists of pion- 
exchange and contact interaction terms. This requires the use of a regularization scheme in order to 
overcome the ultraviolet divergences generated in the momentum integrals when such potentials are 
iterated. The most common approach used to regularize the LS equation is to introduce a sharp or 
smooth regularizing function that suppresses the contributions from the potential matrix elements for 
momenta larger than a cutoff scale, thus eliminating the divergences in the momentum integrals H. 
The second step is to determine the renormalized strengths of the contact interactions, the so called 
low-energy constants (EEC's), by fitting a set of low-energy scattering data. Once the EEC's are fixed, 
the LS can be solved to evaluate other observables. 

Effective field theories and renormalization methods are based on the premise that physics at 
low-energy/long-distance scales is insensitive with respect to the details of the dynamics at high- 
energy/short-distance scales [2 |. In the case of ChEFT, the relevant high-energy effects for describing 
the low-energy observables can be captured in the scale-dependent EEC's of the A^A^ effective inter- 
actions. The A^A^ potential is considered correctly renormalized when the calculated observables are 
approximately independent of the cutoff in the range of validity of the ChEFT. In the language of Wil- 
son's renormalization group, this means that the EEC's must run with the cutoff in such a way that the 
scattering amplitude become (approximately) renormalization group invariant. 

The renormalization group (RG) techniques have been successfully applied to analyze the scale 
dependence and the power counting scheme of A^A^ iteractions in the context of ChEFT |5 6 7| and 
to deiive phase-shift equivalent softer potentials from phenomenological A^A^ interactions by consis- 
tently integrating out high-momentum components, the so called Viowk potentials B8I9I10I . Another RG 
approach that has been recently applied to evolve phenomenological and chiral effective field theory 
(ChEFT) A^A^ interactions to phase-shift equivalent softer forms is the similarity renormalization group 
(SRG) 011I12I13J . The SRG formalism, developed by Glazek and Wilson 014I15II and independently 
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by Wegner lfT6ll . is a renormalization approach based on a series of continuous unitary transformations 
that evolve hamiltonians with a cutoff on energy differences. Viewing the hamiltonian as a matrix 
in a given basis, the similarity transformations suppress off-diagonal matrix elements as the cutoflF is 
lowered, forcing the hamiltonian towards a band-diagonal form and effectively decoupling low-energy 
observables from high-energy degrees of freedom. 

In this work we apply the SRG transformation to evolve an effective A^A^ interaction in leading- 
order (LO) ChEFT, derived within the framework of the subtracted kernel method (SKM) Ill7ll8ll9l20l21l22i 
The SKM is a renormalization scheme in which instead of using a cutoff regularizing function, the LS 
equation is regularized by performing subtractions in the kernel at a given energy scale, while keeping 
the original A^A^ interaction intact. A similar approach based on subtractive renormalization of the LS 
equation is described in Ref. 1.23.24,1 . although a momentum cutoff is also introduced to regularize the 
momentum integrals. 



2 Subtractive Renormalization 

In the following, we show how to derive the subtracted kernel LS equation ifTTl . Here and in what 
follows we use units such that Ti - c - M - 1, where M is the nucleon mass. 
Consider the formal LS equation for the T-matrix, written in operator form as 

T{E) = y + y GliE) T(E) = y [l + G,;(£) T(E)] , (1) 

where V is the interaction potential and Gg (£) is the free Green's function for the two-nucleon system, 
given in terms of the free hamiltonian Hq by 

G;(E) ^ [E - Ho + ie]-' . (2) 

Using Eq. ([T]i, we can express the potential V in terms of the T-matrix at a given energy scale -/u^, 

V = r( V) [l + G(-f/)T(-f/)l' . (3) 
Putting Eq. ^ back in Eq. ([T]i, we obtain 

+ T(E)G;(E)T(-i/) [l + G;(-i/)T(-f?)l' . (4) 

Now, we can multiply the entire equation by the term in the square brackets to find the subtracted 
kernel LS equation, 

T(E) = T(-fi^) + [g;(E) - G,;(-//2)j T(E) , (5) 

G<'>(£;-/j2) 

The Green's function after one subtraction, G^^\E; -ju^), can be written as 

G^'\E; -f?) ^ \{-^? - E) G,t( V)! G^iE) = I'^'f). Gl{E). (6) 

The subtracted kernel equation has the same structure as the formal LS equation Eq. ([T]i but with the 
potential y replaced by Ti-fi^) and the free Green's function, Gq{E), replaced by G^^\E; -//^), which 
contains one subtraction at the energy scale -ju^. For a pure contact interaction, the subtracted kernel 
LS equation produces a finite T-matrix due the presence of the form-factor in G^^K 
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2.1 Recursiveness and Multiple Subtractions 

If we consider an interaction that contains stronger divergencies, like the ChEFT A^A^ interaction up 
to NNLO, we need more subtractions to render the T-matrix finite. As the degree of the divergency 
increases, the kernel subtractions can be performed recursively to obtain a finite T-matrix I.18J. The 
general LS equation for multiple recursive subtractions is given by 

Tn{E) = V^"\E; -I?) + y(">(£; -i/)G'^"\E; -i/)TiE) , (7) 

where the generic n-subtracted kernel is 

G^"\E; -^?) = [{-t? - E) Gti-t?)]" Gl(E) = ^j^^f^ Gt(E). (8) 

The recursive driving terms can be written as 

y<"'(£; -i?) = y<"-'^(£; -^?) + y<"-"(£; -^?) g'"\E; -^^) y^"' , (9) 

with 

_^2) ^ (_^2 _ [g,;(-^2)1" . ^=4^^ . (10) 

i-f^ - no)" 

In order to obtain the n-subtracted T-matrix from Eq. (|7]i, first one needs to recursively solve Eq. (|9]l 
for y*^"' up to the number of subtractions required to regularize the considered interaction. 



2.2 Fixed-point Interactions and Renormalization Group Invariance 

Given the driving term V^"\ we can construct a fixed-point hamiltonian lfT9]| . = Ho + V-r, where Vk 
is the renormalized interaction. Both /f^ and are fixed-point operators, i.e. invariant with respect to 
the subtraction scale -fp-. 

Replacing Vk in Eq. ([T]i we obtain the LS equation for the corresponding renormalized T-matrix, 

TniE)^VK + VkG(E)Tk(E). (11) 

By construction, the T-matrix obtained from the solution of this equation must be equivalent to that 
obtained from the solution of the n-subtracted kernel LS equation: T^iE) = T„(E). From this condition 
we obtain an integral equation which relates the renormalized interaction Vk to the driving term V^"\ 

= V^"\E; -f/) - y(«>(£; -^u^) _^2^ ^ (12) 

with 

(V-^r 



(n) 
9k 



1 



(-Ai2 - Ho)" 



G^o(E) ■ (13) 



The subtraction scale -fr is arbitrary and so all observables should not depend on its choice. This 
condition is fulfilled by imposing the invariance of the T-matrix with respect to -fi^, which yields to a 
renormalization group equation in the form of a non-relativistic Callan-Symanzik (NRCS) equation. 



V"\e- V) = -yw "\ ; ^ ^ y<")(£; -ju2) . (14) 



with the boundary condition V^"\E; -yL(^)|^^/i - V^"\E; -fP') imposed at some reference scale jl 
where the physical information is supplied. This equation governs the evolution of the driving term 
V^"\E; -fp-) with the subtraction scale -f/ in such a way that the T-matrix remains invariant. 
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Fig. 1. Left panel: Evolution of V^^\k, k) with the subtraction scale -yu^; Right panel: Relative differences between 
the phase-shifts in the '5o channel calculated at a scale yu = 10 fm"' and at yu — > oo in three distinct cases. 



In order to illustrate the renormalization group invariance in the SKM approach, we consider the 
LO ChEFT A^A^ interaction in the 'S'o channel, which consists of the one-pion exchange potential 
(OPEP) plus a Dirac-delta contact interaction. The driving term for the corresponding one-subtracted 
LS equation (n = 1) is given by 



y('>(£;V) = VoPEP + Co(V) 



(15) 



where Co(-/U^) is the renormalized strength of the contact interaction. 

We solve Eq. (14 1 for the matrix elements y*''(/?,/?') of the driving term projected in the '5o 
channel, using a partial-wave relative momentum space basis. The boundary condition is set at the 
reference scale /i — > oo, where the renormalized strength of the contact interaction Co(-jU^) is fixed 
by fitting the experimental value of the scattering length in the '^o channel, a, - -7^1 fm. Once the 
strength Co(-A(^) is fixed, and so the matrix elements y*'^(p,p') of the driving term are known, we 
can compute the coiTesponding matrix elements of the renormalized interaction Y^ip, p') by solving 
Eq. (12 1 numerically. In the left panel of Fig. [T] we display the result obtained for the evolution of the 
diagonal matrix element V'-^\k,k) with the subtraction scale -p^ for k = - 1-55 fm"'. As one 
can observe, V^^\k,k) is enhanced (becoming more attractive) as p decreases. In the limit — > oo, the 
driving term becomes independent of p and matches the corresponding renormalized potential. 

Renormalization group invariance can be easily verified by evaluating the phase-shifts in the '5o 
channel as a function of the laboratory energy ^lab- In the right panel of Fig. ([T]l we show the relative 
differences between the phase-shifts calculated at p - 10 fm"' and at ju — > oo in three distinct cases. 
The black line corresponds to the result obtained from the solution of the subtracted kernel LS equation 
with the driving term V^'* determined by simply fixing the strength of the contact interaction Coi-p^) 
at p - 10 fm"' to fit the scattering length. The red line corresponds to the result obtained by evolving 
the driving term V''* through the NRCS equation from p — > ootoyu - 10 fm"'. The blue squares 
correspond to the result obtained by solving the LS equation for with the renormalized potential 
V-R determined from the driving term V*'^ evolved through the NRCS equation. As one can observe, 
when the phase-shifts are evaluated with the unevolved driving term there is a residual dependence on 
the scale p due to the fitting procedure used to fix the renormalized strength Coi-p^). On the other 
hand, when the driving term is evolved through the NRCS equation the phase-shifts remain invariant, 
apart from relative numerical errors smaller than 10"^. One should note that the phase-shifts obtained 
from the LS equation for the renormalized T-matrix are also invariant, since the renormalized potential 
determined from the evolved driving term is a fixed-point operator. 
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3 Similarity Renormalization Group for A^A^ Interactions 
3.1 Theoretical background 

The general formulation of the SRG approach was developed by Glazek and Wilson [14' 151 in the 
context of light-front hamiltonian field theory, aiming to obtain effective hamiltonians in which the 
couplings between high and low-energy states are eliminated, while avoiding artificial divergences 
due to small energy denominators. 

Consider an initial hamiltonian in the center of mass frame for a system of two nucleons, which can 
be written in the form H - Ty^i + V, where r,ei is the relative kinetic energy and V is the A^A^ potential. 
The similarity transformation is defined by a unitary operator designed to act on the hamiltonian and 
evolve it with a cutoff A on free energy differences at the interaction vertices, 

H, = UiA) H u\A) = r,ei + /, Va , (16) 

where /i is a similarity function and Va is called the reduced interaction. The unitarity condition 
satisfied by U(A) is given by: 

U(A) uUA) = U\A) U(A) = 1 . (17) 



The similarity function /a is a regularizing function which suppresses the matrix elements between 
states with free energy differences larger then the cutoff A, such that the hamiltonian is driven towards 
a band-diagonal form as A is lowered. Usually, /a is chosen to be a smooth function of the similarity 
cutoff A. A simpler choice is to use a step function. The similarity transformation can be defined in 
terms of an anti-hermitian operator jja which generates infinitesimal changes of the cutoff A, 

m - UiA) ^ = -.1 . (18) 

Using this definition and the unitarity of U{A) we can derive a first-order differential equation for the 
evolution of the hamiltonian, 

^ = (19) 

with the boundary condition = H. 

In the application of the SRG described in this work, we employ the formulation developed by 
Wegner II61I . based on a flow equation that governs the unitary evolution of the hamiltonian 

^^[ns,H,], (20) 
as 

with the boundary condition //jLq = H. The flow parameter s has dimensions of (energy)"^ and 
ranges from to 00. In terms of a similarity cutoff A, here with dimensions of momentum, the flow 
parameter is given by the relation s = A^'^. 



Wegner's flow equation is analogous to Eq.( 19 1, but the specific form rj^ - [Gj, i/j] is chosen for 



the anti-hermitian operator that generates the unitary transformation, which gives 

^ = [[G„HA,H,] . (21) 
as 

Such a choice for the generator rj, corresponds in the Glazek- Wilson formulation to the choice of a 
gaussian similarity function /) with uniform width A . The operator Gs defines the generator rjs and so 
specifies the flow of the hamiltonian. Wegner's choice in the original formulation is the full diagonal 
part of the hamiltonian in a given basis, G , = diagiH,). A simpler choice is to use the free hamiltonian, 
Gj = Trei, which gives the generator tJs - [Ti-ei, Hs] 
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Fig. 2. (Color online) SRG evolution of the SKM-LO ChEFT potential in the '5o channel (in units of fm). 



3.2 SRG Evolution in the SKM Approach 

Using the generator t], - [Tre\, H^], Wegner's flow equation for the SRG evolution of the A^A^ potential 
matrix elements is given by 

f'^ ^ -ip^-p'^fV,ip,p') + - dqqUp^+p'^-2q^)VAp,q)Vsiq,p'). (22) 
as n Jo 

For simplicity, we are using Vs{p, p')asa short notation for the projected A^A^ potential matrix elements 
y(jLLS,i)^^^ p') in a partial-wave relative momentum space basis, | q{LS )J; I ), with normalization 



2 r 

1 = - dqq^\q(LS)J;I}{q(LS)J;I\, (23) 
^ Jo 

where the indexes J, L{L'), S and / respectively denote the total angular momentum, the orbital angular 
momentum, the spin and the isospin quantum numbers of the A^A^ state. For non-coupled channels 
{L - L' - J), such as the singlet '^o, the A^A^ potential matrix elements Vs(p,p') are simply given by 
Vs(p,p') - V\^''^^ ''\p, p'). For coupled channels {L,L' = 7 ± 1), such as the triplet D\, the 

Vsip, p') represent 2x2 matrices of matrix elements for the diflTerent combinations of L and L', 

y(pp') = (^^""''^^'^'^^^""'"^^'^'^^ (24) 

^ sKP^P ) - \ JJL'LSJ), ,^ ^AJL'L'S;!), ,J- ^^^^ 

As a consequence of the choice of the SRG transformation generator rjs - [T^ei,Hs] each interaction 
channel evolves with the cutoff' A - 5 '^"^ independently of the other channels |25J. 



We solve Eq. (22i numerically, obtaining an exact (non-perturbative) solution for the evolution 
of the SKM-LO ChEFT potential in the '^o channel. The relative momentum space is discretized 
on a grid of A^ gaussian integration points (we have used 200 mesh points), leading to a system of 
A^^ non-linear first-order coupled differential equations which is solved using an adaptative fifth-order 
Runge-Kutta algorithm. The boundary condition is set at s = (/I — » oo), such that the initial potential 
is given by the fixed-point renormalized interaction V'ji(p,p') derived through the SKM scheme. As 
one can observe from Fig. (j2]l, the off-diagonal matrix elements are systematically suppressed as the 
similarity cutoff A is lowered, such that the potential is driven towards a band-diagonal form. 
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Fig. 3. Top panel: SRG evolution for the renormalized potential; Bottom panel: Evolution of the driving term 
through the NRCS equation. Units are fm for the potential and frrT^ for the momenta. 



In Fig.[3]we compare the SRG evolution of the matrix elements Vk(p,p') for the fixed-point renor- 
malized interaction with the similarity cutoff A (top) and the evolution through the NRCS equation of 
the matrix elements y'^'*(p, p') for the driving term with the subtraction scale yu (bottom). As one can 
observe, the NRCS evolution preserves the global shape of the driving term. On the other hand, the 
SRG evolution completely modifies the shape of the renormalized interaction. 

In the left panel of Fig. (|4| we show the phase-shifts in the '5o channel obtained for the initial 
fixed-point renormalized interaction V-Rip, p') and for the corresponding SRG potentials evolved up to 
several values of the similarity cutoff A. As expected for a unitary transformation, the results are the 
same (apart from relative numerical errors smaller than 10"^). 

We now test for the decoupling of low-energy observables from high-energy degrees of freedom 
following the analysis introduced by Bogner et al. 11121131 . which consists in applying an exponential 
function to the SRG potential that suppresses contributions from matrix elements Vs{p,p') with p,p' 
larger than a given momentum k^^^, 

Vf—"\p, p') = exp[-(p2/^^^J"] VAp, p') exp[-(p'2/fc2 ^,)«] , (25) 

In the right panel of Fig. Q we show the phase-shifts obtained by cutting the initial potential V'Rip, p') 
and the SRG potential evolved to a similarity cutoff A - 2 fm"' at fcmax = 3.5 fm"', with n = 8. As one 
can observe, the phase-shifts obtained for the cut initial potential V<r{p, p') are completely modified 
in comparison to those for the uncut potential. For the cut SRG potential, the phase-shifts remain 
practically unchanged at low energies. Therefore, an explicit decoupling between the low- and high- 
momentum components is verified for the SRG evolved potential. 



4 Summary and Concluding Remarks 

We have investigated the similarity renormalization group evolution (SRG) of A^A^ interactions in the 
framework of the subtracted kernel method (SKM), a renormalization scheme based on subtractions 
performed in kernel of the scattering equation. 
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Fig. 4. Left panel: Phase-shifts in the '5o channel as a function of the laboratory energy; Right panel: Test for the 
decoupling between low- and high-momentum components. 



Considering a simple example, the LO ChEFT A^A^ interaction in the '5o channel, we have shown 
that a fixed-point renormalized interaction and renormalization group invariant phase-shifts can be 
obtained from the subtracted kernel scattering equation, provided the driving term is evolved with the 
subtraction scale through a renormalization group equation in the form of a non-relativistic Callan- 
Symanzik (NRCS) equation. We have solved Wegner's flow equation numerically to obtain a non- 
perturbative solution for the SRG evolution of the fixed-point renormalized interaction. By calculating 
the phase-shifts, we have verified the unitarity of the similarity transformation. By cutting the SRG 
potential at a given momentum using an exponential function, we have verified the decoupling of 
low-energy observables from the hig-energy degrees of freedom. 

The next step, to be implemented in a future work, is to consider the the SRG evolution of A^A^ 
interactions up to higher-orders in ChEFT and other partial-wave channels. 
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